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ABSTRACT 


To study radiation effects caused Lv geomagnetically trapped high energy 
protons, one must produce a proton source with the expected spectral charac- 
teristics. A method for generating such spectra using a spinning wheel of 
varying thickness is described. Resulting wheel shapes for two typical proton 
spectra are shown. 
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A METHOD FOR SIMULATING VAN ALLEN BELT 
PROTON ENERGY SPECTRA 

INTRODUCTION 


Film stacks will be carried to observe cosmic ray interactions on 
certain space flights passing through the magnetically trapped radiation belts 
around the earth. Although, in general, these flights will remain in the belts 
only short periods, some determination is needed of the expected film degrada- 
tion caused by trapped protons. In making this determination, it is necessary 
to simulate the trapped proton energy spectrum by using a proton accelerator 
and irradiate the film to be used. It is extremely difficult, if not impossible, 
to vary the energy or number of protons rapidly in a controlled fashion for *aiy 
of the accelerators now in use. The problem, then, is to establish a 'iethod 
that will allow rapid variation in the proton number and the energy of a proton 
accelerator beam; the goal is to simulate a given proton spectrum that can vary 
in number over an order of magnitude in the energy range between 10 MeV and 
200 MeV. 


THE METHOD 


In the energy range of interest, two approximations (the straight-ahead 
approximation and the continuous slowing-down approximation) are valid for 
describing the transport of protons through material shields. Stated simply, 
these approximations mean that a proton will not be scattered from its initial 
direction by the shield and that the energy of the proton at any point in the shield 
is a continuous function of the initial energy and the distance traveled through 
the shield. Thus, if a slab of material is placed in a proton beam, the beam 
will come out of the other side, traveling in the same direction and containing 
the same number of particles, but the beam will have a reduced energy that 
will depend upon the initial beam energy, the thickness of the slab, and the type 
of material in the slab. By placing the film in such a beam described above and 
varying the slab thickness correctly with time, the effects of exposing the film 
to a given spectrum can be simulated. One process to quickly vary the slab 
thickness is to use a spinning wheel of varying material thickness as shown in 
Figure 1. 




FIGURE 1. THE PROTON SPECTRUM WHEE L CONFIGURATION 


ASSOCIATED PROBLEMS 


The majority of the energy lost by a proton in passing through a slab is 
through interactions with the electrons in the slab. This is a statistical process; 
consequently, the protons coming out of the slab would be expected to have some 
distribution of energies instead of having a discrete energy. This corresponds 
to a distribution of path lengths about the mean path length. Reference 1 
tabulates the root mean square straggling (the standard deviation of the distri- 
bution; for the mean path length required for a proton to lose all its energy. 

(This path ie.igth is the length of the actual path followed as opposed to the range, 
which is the straight line distance traversed in stopping. ) For a 200 MeV proton 
stopped in aluminum, the deviation is 0. 36613 g/cm 2 about a mean path length of 
33. 201 g/cm 2 . For the range, which is 33. 112 g/cm 2 , the deviation should be 
about the same since the range is so close to the mean path length. The range 
of a 15-MeV proton is 0. 34297 g/cm 2 . Thus, if a slab of range thickness is 
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placed in front of a 200 MeV proton beam, about half of the protons will be 
stopped, and most of those getting through the slab will have energies less than 
15 MeV. This is the worst case possible. For thinner slabs the deviation of 
the energy distribution will be considerably less than 15 MeV. Since the trapped 
belt proton spectra to be reproduced do not have any sharp peaks, the smearing 
effects caused by the inability to produce discrete energies are not as important 
as they seem at first. 

Perhaps the main problem in using a materj^l slab to reduce a proton 
beam’s energy is the possibility of nuclear interactions in the slab; this would 
result in the scattering of protons out of the beam, production of neutrons, and 
activation of the slab. The use of some light material such as aluminum will 
prever t activation of the slab. 

The nuclear interactions can be divided into two classes, elastic and 
inelastic collisions. Elastic collisions leave the nucleus unexcited and produce 
no secondary particles, but the proton is scattered at some angle with its 
original direction. Generally, this scattering is very peaked in the forward 
direction except for materials of very low atomic number. In inelastic collisions 
for the energy range under consideration, the proton loses part or all of its 
energy to the nucleus, one or more secondary protons and neutrons may be 
produced, and again the proton is lost from the beam. R. G. Alsmiller of 
Oak Ridge National Laboratory, using a high-energy nucleon Monte Carlo 
transport code, calculated the expected neutron and proton spectra that would 
be produced by 160 MeV protons incident on an aluminum slab [2]. Some of 
the results of this calculation are shown in Table I. As can be seen, secondary 
protons can never account for more than about 3 percent of the total flux, and 
secondary neutrons make up about 20 percent of the total. The fractions given 
are for the total flux and not just the flux in the dir ction of the beam. A large 
fraction of these secondaries are expected to be scattered out of the beam 
direction or emitted isotropically in the case of evaporation neutrons. 

In Table I the secondary neutrons are divided according to the process 
that produced them. The cascade neutrons, which account for about 5 percent 
of the total flux, are produced in direct interactions with the incident proton 
inside the nucleus. They are fairly energetic and have a forward-peaked 
angular distribution. Evaporation neutrons, which account for about 15 percent 
of the total flux, are produced by a process in which a proton excites the target 
nucleus, which, in turn, loses its excitation energy by the emission of a neutron 
or proton. These particles have fairly low energies ( <20 MeV) , and their 
angular distribution is isotropic. The angular distribution of evaporation 
neutrons causes that component of the flux to decrease very rapidly with increas- 
ing distance from the slab. (Most of the evaporation protons never get out of 
the slab. ) 
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TABLE I. PRIMARY AND SECONDARY PROTON TRANSMISSION 
PER UNIT FLUX THROUGH VARIOUS ALUMINUM SHIELD THICKNESSES 


Shield 

Thickness 

(g/cm 2 ) 

Primary 

Protons 

Secondary 

Protons 

Cascade 

Secondary 

Neutrons 

Evaporation 

Secondary 

Neutrons 

5 

0.9525 

0.018 

0.028 

0. 100 

1° 

0.9025 

0.030 

0.048 

0. 126 

15 

0.8585 

0.024 

0.063 

0. 170 

21 

0.7975 

0.001 

0.068 

0. 180 


The mean free path of a 20-MeV neutron between nonelastic collisions is 
18 cm of aluminum. Thus, very few neutrons would be expected to have an inter- 
action in a piece of photographic film. They might produce some activation in 
the surrounding equipment, but the fluxes are so small that this should not be a 
problem. If the total flux is not accounted for in the simulation, the secondary 
protons would be expected to slightly raise the high-energy side of the spectrum 
produced, since the low energy secondary protons are stopped in the slab. It 
was concluded that this was net a significant problem. By knowing the cross 
section for nuclear interactions for protons, it is fairly easy to determine the 
fraction of the primary protons scattered out of the beam and correct for this in 
calculating the thickness variation around the circumference of the wheel. 


METHOD FOR CALCULATION OF THE WHEEL THICKNESS VARIATION 


To determine the variation in the wheel thickness to produce a given spec 
trum given an incident beam of protons with energy and flux I Q , we proceed as 
follows. If the differential proton spectrum is described by 0(E), then 


df = 


0 (E) dE 



( 1 ) 


where df is the fraction of the total proton flux between 0 and E 0 that has energies 
between E and E + dE. If dS (E) is the fraction of the wheel circumference that 
is thick enough to reduce E 0 to energy E, then 
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( 2 ) 


df = 


I n c~ 2T<Eo ’ E) dS(E) 


f 1 . - ST(E„ E) 

J i 0 c 


d S 


where T(E 0 , E) is ihe wheel thickness required to produce energy E from E 0 ; 

—i- T ( E E) 

the exponential factor e ~ i() ' is to correct for protons lost from the 

beam by nuclear interactions; Z is the nuclear inelastic cross section; and the 
integral is just the total primary flux coming out of the slab. From equations 
( 1) and (2), 


0 (E) dE 


, E . 

J 4 (E) d E 


c- ETlE - E| dS(E. 


S' i, „- rTIE » E 'd, 


(3) 


or 


^.- rT,E>E> dS 

( 0(E) dE 


It is known that 


/' d S (E) = 1 ; 


(5) 


thus, 


e - £T(E »' E| dS A(E)o £T(E »' E) dE 


J 0(E) dE 


= 1 ( 6 ) 


or 


I 


I 


1 e - s T (E 0 , E) d s = 0 (E) dE 


J 0 (E) e 
J o 


2T (E 0 , E) 


d E 


(7) 


From equations (4) and (7), one gets 


5 



dS (E) 


( 8 ) 


0 (E) C LT(E °’ E) d E 
/ E ° 0 <E)e 2;T < E - E) dE 


cr 


S(E’) = 


l 0 (E)e ST(E »’ E) dE 
J E ° 0 (E) e 2T(E «* E) dE 


(9) 


To obtain T(E 0 , E) , one uses the relationship 

R(E 0 ) = R(E) + T(E 0 , E) (10) 

or 

T (E 0 , E) = R (E 0 ) - R (E) , (11) 

where R(E) is the mean range of a proton at energy E. 

To minimize problems caused by secondaries and activation, aluminum 
was selected as the material for the wheel, and a 160-MeV beam energy was 
used. The proton range was approximated by 

R(E) = ^ in (1+ 2bE r ) (12) 

from Reference 3. /or aluminum, values for the parameters a, b, and r are 

—3 —6 

as follows: a = 2.794 x 10 , b = 2. 346 x 10 , and r = 1. 775 from Reference 4, 

where E in in millions of electron volts and R(E) is in grams per square centi- 
meter. (The fit is good to better than 3 percent ove*' the range of interest. ) 

The aluminum density assumed was 2.692 g/cm 3 . 

Since the inelastic cross section is a funcvion of energy and since the 
energy of the proton changes as it passes through Jie slab, the apparent inelas- 
tic cross section to account for the primary protons lost in a slab is a function 
of the initial proton energy and the slab thickness. To obtain the effective in- 
elastic cross section, Alsmiiler's Monte Cano results [2] for primary flux 
transmission were taken from Table 1, the effective cross section was deter- 
mined, and these values were combined with the known cro3S section for 160-MeV 
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protons and fitted with the equation. 

-4 -3 

2 = 4.01 x 10 T+ 8.92 x 10 , 

where T is the thickness in centimeters, and 2 is the effective inelastic cross 
section in square centimeters per gram. Elastic scattering was assumed to be 
so peaked in the forward direction that in these calculations it could be ignored. 

RESULTS 

Figure 2 shows the wheel thickness, V, in centimeters versus the frac- 
tion of the total circumference, S, having thickness greater than T. The spec- 
trum to be simulated, shown in Figure 3, was obtained using data described in 
Reference 5 and a computer program described in Reference 6. 



FIGURE 2. WHEEL SHAPE FOR THE SPECTRUM OF PROTONS IN A 
389 km (210 n. mi.), 50 DEGREE CIRCULAR ORBIT 
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The actual wheel used in the experiment was slightly more complicated 
than the one shown in Figure 2. The incident beam had an energy of 160 MeV, 
but there were plastic scintillators 1. 296 g/cm 2 thick located in front of and 
behind the wheel to monitor the proton beam. Thus, the highest energy passing 
through the whole apparatus was 145.9 MeV. The proton spectrum to be simul- 
ated was known to 600 MeV. To account for particles with energies greater than 
145.9 and less than 600 MeV, two gaps were left in the wheel. Their total width 
was proportional to the fraction of the total flux with energies between 145.9 MeV 
and 600 MeV. The resulting wheel shape is shown in Figure 4. The spectrum 
to be simulated is shown in Figure 5. 
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FIGURE 4. WHEEL SHAPE FOR THE SPACTRUM SHOWN IN FIGURE 5 
FOR AN INCIDENT BEAM OF ENERGY OF 160 MeV, 

WHERE PLASTIC SCINTILLATORS 1.296 g/cm J 
ARE PLACED BEFORE AND AFTER THE WHEEL 






Lester Katz did the actual mechanical design work on the machining of 
the wheel shape and the development of a drive mechanism to spin the wheel. 
Figure 6 shows a photograph of the final operational device used in the exper- 
iment. 



FIGURE 6. THE WORKING APPARATUS USED IN THE EXPERIMENT 
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